arXiv:l502.01832v2 [astro-ph.GA] 7 Mar 2015 


Mon. Not. R. Astron. Soc. OOO.fTp?? (2011) Printed March 10, 2015 (MN [AT E X style file v2.2) 


Neural-network selection of high-redshift radio quasars, and the 
luminosity function at z ~ 4 


D. Tuccillo, 1 * J. I. Gonzalez-Serrano , 1 C. R. Benn. 2 

1 Instituto de Fisica de Cantabria (CSIC - Universidad de Cantabria), E-39005, Av. de los Castros s/n Santander, Spain 
2 Isaac Newton Group, Apartado 321, E-38700, Santa Cruz de la Palma, Spain 


accepted for publication in MNRAS on 3 March 2015 


ABSTRACT 

We obtain a sample of 87 radio-loud QSOs in the redshift range 3.6 ^ z ^ 4.4 by cross- 
correlating sources in the FIRST radio survey (Siaghz > 1 m Jy) with star-like objects having 
r < 20.2 in SDSS Data Release 7. Of these 87 QSOs, 80 are spectroscopically classified 
in previous work (mainly SDSS), and form the training set for a search for additional such 
sources. We apply our selection to 2,916 FIRST-DR7 pairs and find 15 likely candidates. 
Seven of these are confirmed as high-redshift quasars, bringing the total to 87. The candidates 
were selected using a neural-network, which yields 97% completeness (fraction of actual 
high-z QSOs selected as such) and an efficiency (fraction of candidates which are high-z 
QSOs) in the range of 47 to 60%. We use this sample to estimate the binned optical luminosity 
function of radio-loud QSOs at 2 ~ 4, and also the LF of the total QSO population and its 
comoving density. Our results suggest that the radio-loud fraction (RLF) at high z is similar 
to that at low-z and that other authors may be underestimating the fraction at high-z. Finally, 
we determine the slope of the optical luminosity function and obtain results consistent with 
previous studies of radio-loud QSOs and of the whole population of QSOs. The evolution 
of the luminosity function with redshift was for many years interpreted as a flattening of 
the bright end slope, but has recently been re-interpreted as strong evolution of the break 
luminosity for high-z QSOs, and our results, for the radio-loud population, are consistent 
with this. 
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1 INTRODUCTION 


Since their discovery in the 1960s ( [Hazard, Mackey & Shimmins| 
|1963||Schmidt|l963) , quasi-stellar objects (hereafter QSOs) have 
played a key role in extragalactic research, in particular in connec¬ 
tion with studies of super-massive black holes (hereafter SMBHs), 
galaxy evolution, the intergalactic medium, large-scale structure 
and cosmology. 

Quasar candidates are mainly identified from their optical 
colours in large sky surveys such as the 2dF survey ((Boyle et al.| 
20001) and the Sloan Digital Sky Survey (SDSS, |York et al.|2000[ 


Schneider et al7| |20 1 0[ . Current data from SDSS provide us with 
photometric measurements for « 5 x 10 8 galaxies, quasars and 
stars. The survey also provides spectra for nearly two million 
of these objects. This has dramatically increased the number of 
known QSOs, since the first edition of the SDSS quasar catalogue 
( jSchneider et al.|20Q2( l. The 5th edition, which is based upon the 
SDSS Seventh Data Release(DR7), includes a remarkable 105,783 
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spectroscopically-confirmed QSOs ( jSchneider et al.||~2010[ l. Most 
of the SDSS QSO candidates were selected as spectroscopic tar¬ 
gets on the basis of their non-stellar colours in ugriz bands and 
by matching unresolved sources to the FIRST radio survey (Faint 
Images of the Radio Sky at Twenty-Centimeters,(Becker, White &| 
|Helfand|1995) . 

The SDSS QSO-selection algorithm was presented by 
Richards et al. (2002), and according to them is sensitive to QSOs 
at redshifts z <~ 5.8. Completeness (fraction of QSOs selected 
as such) and efficiency (number of actual QSOs amongst the candi¬ 
dates, divided by the total number of candidates) of the selection are 
a complex function of apparent magnitude i and redshift. Although 
QSOs of type 2 and certain QSOs of type 1 are missed, the overall 
estimated completeness is high (Van demBeAet al.|2005) , above 
90% for 16.0 < i < 19.0 ( [Richards et al.|2006( . At higher red¬ 
shift, both completeness and efficiency drop, with an overall com¬ 
pleteness of ~ 80% for 3 ^ z ^ 5.3 and efficiency ~ 55% for 
QSOs with 2 > 3. 

Determination of the QSO Luminosity Function (QLF) is im¬ 
portant for the study of active galactic nuclei (hereafter AGN) and 



























2 D. Tuccillo et al. 


it requires QSO samples with a statistically-significant number of 
sources and with accurately-known completeness. In particular an 
accurate knowledge of the QLF at different epochs allows impor¬ 
tant constraints to be placed on the evolution of the mass function 
of SMBHs, on their growth and on the lifetime of the QSO phase 


i jVolonteri, Haardt & Madau]|20031 Netzer & Trakhtenbrot||2007| 


; Shankar et al.||20l0 ; Shen & Kel y 2012}. The QLF also pro 
vides important information about the impact of QSO activity on 
the formation and evolution of the host galaxies |Cattaneo et al. 
|2009| [Fabian 2012). In addition it allows constraint of the contri¬ 
bution of AGN to the X-ray background (|Ueda et al.|2003||Hickox 


& Markevitch|2006}, the ultraviolet ionizing radiation )|Samantaray 


&KhareJ2000| Giallongo et al.|2012| > and the infrared background 
( |Dole et al.|2006) . 

Radio-loud QSOs (hereafter RLQs) account for ~ 8 — 13% 
of QSOs (Ivezic et al.||2002| |Jiang et al,||2007| |Balokovic et al.| 
|2012[ ). The exact radio-loud fraction (RLF) is still unknown and 
it also depends on the definition of radio-loudness. Usually the 
latter is based on the radio luminosity P emitted by the source 
(e.g. Gregg et al. | |1996} , define an object to be radio-loud if 
log Pi.4,GHz.(W/ Hz) > 25.5), or on R , the ratio between 

monochromatic radio and optical luminosities ) |Stocke et a1. 1 1992} . 
Some authors see no evidence for significant change of RLF with 
either redshift or luminosity (e.g., Goldschmidt et al. 1999] |Stem| 


|et al.|2000}|Cirasuolo et al.|2003[|Vigotti et al.|2003 1 , while others 

find that the RLF changes with both parameters (e.g.|Miller, Pea^| 


cock & Mead|1990[ Visnovsky et al.|1992[|Schneider et al.|1992| 


Jiang et al.|2007||Balokovic et al.|2012). In particular! Jiang et al.| 

| 2007} 7 using a sample of 30,000 optically selected quasars from 
the SDSS quasar catalogue obtained from Data Release 3 l |Schnei-| 
der et al. 2005), find that the RLF of quasars decreases strongly 
with increasing redshift and decreasing luminosity. 

Although only a small fraction of QSOs are radio-loud and, 
while they may not be fully representative of the entire population, 
they can be used to check the completeness of samples selected 
mainly on the basis of optical data, such as the SDSS QSO cata¬ 
logues. For example, selection at radio wavelengths greatly reduces 
incompleteness due to the effects of dust obscuration, reddening 
and/or the presence of broad absorption line s fCarballo et al.|2008| 
hereafter C08, and |McGreer, Helfand & White|2009) . Selection in 
the radio is also more efficient, since it reduces stellar contamina¬ 
tion. 

In addition, increasing the number of known radio-loud QSOs 
is important for understanding the origin of the radio phenomenon 
itself and for clarifying the connection between radio and optical 
activity in QSOs. Indeed, whether RLQ and radio-quiet QSOs are 
two physically distinct populations (the so-called radio-loudness 
dichotomy) or are part of a continuous sequence is still a matter 
of debate (see |Jiang et al.||2007[ |Balokovic et al. 2012| |Kratzer| 
2014|>. As pointed out by several authors (e.g. Cirasuolo et al.|2003[ 


Balokovlc et al.|2012|>, it would be helpful in the context of this de¬ 


bate, to have more-homogeneous samples of QSOs, a larger num¬ 
ber of QSOs with radio data available, and a reduction of selection 
biases in flux-limited samples. 

In this work we present the results of a selection of RL QSOs 
with 3.6 ^ z ^ 4.4. We start by cross-matching FIRST radio 
sources (from the 2003 April 11 release of the catalogue), with star- 
like objects in the SDSS DR7 photometric catalogue (Section 2). 
To select candidate high-redshift QSOs (Section 3), we use a sim¬ 
ple but reliable neural-network (NN) algorithm already tested in 
previous work jCarballo et al.|2006||Carballo et al.|2008} . Spectra 
of the resulting 15 candidates were obtained with the 2.5m Nordic 


Optical Telescope (NOT) on La Palma (Spain) and these are pre¬ 
sented in Section 4. together with checks of the completeness of 
our selection and comparisons of the completeness and efficiency 
with other work. In Section 5 we discuss the K-correction, and we 
present our sample of 87 3.6 ^ z ^ 4.4 radio-loud QSOs. The 
various sources of incompleteness of this sample are discussed in 
Section 6. In Section 7, we calculate the optical luminosity func¬ 
tion for 2 ~ 3.8 and 2 ~ 4.2. In Section 8 we derive the space 
density of radio-loud QSOs, and the space density of all QSOs, and 
discuss the luminosity function. Our conclusions are summarized 
in Section 9. 

All optical magnitudes are on the AB system. We use a 
ACDM cosmology with Da = 0.7, = 0.3, and Ho = 

70 km s -1 Mpc -1 . 


2 DATA 

2.1 Surveys used 

In this work we use the FIRST radio survey and the SDSS DR7 
optical survey to obtain a radio-optical sample of QSO candidates 
in the redshift range 3.6 ^ 2 ^ 4.4 . 

The FIRST survey was originally designed to produce the ra¬ 
dio equivalent of the optical Palomar Observatory Sky Survey, us¬ 
ing the NRAO Very Large Array (VLA) in its B-configuration at 
1.4 GHz. Subsequently the survey area was chosen to make it ideal 
for comparison with the Sloan Digital Sky Survey (SDSS). The sur¬ 
vey produces images with 1"8 pixels, a resolution of 5", a typical 
rms of 0.15 mJy and a flux-density limit of 1 mjy. The positional 
accuracy at the survey flux limit is ~ 1". We used the 2003 April 11 
version of the FIRST catalogue containing 811,117 sources cover¬ 
ing a total area of 9,033 deg 2 (8,422 deg 2 in the Northern Galactic 
Cap and 611 deg 2 in the Southern Galactic Cap). 

The SDSS DR7 JAbazajian et al.|2009} covers a total imaging 
area of 11,663 deg 2 (7,646 deg 1 in the Northern Galactic Cap). 
A total of 357 million distinct objects are included in the imaging 
catalogue, of which approximately 1.6 million are also included in 
the spectroscopic catalogue. 

The survey reaches magnitude limits (95% detection repeata¬ 
bility for point sources) in photometric bands u, g, r, i and 2 of 
22.0,22.2,22.2,21.3 and 20.5 respectively. Absolute astrometric 
errors are < 0 , /l. In this paper, we consider only the images flagged 
by SDSS as ’Primary’. These are unique detections, i.e. they do 
not include duplicate detections from the overlap between survey 
stripes. Each such object is associated with a run and a field which 
is the primary source of imaging data at this position. 

To determine the overlap area of the two surveys we first deter¬ 
mined the area of the FIRST survey, which has an irregular bound¬ 
ary, by constructing a Delaunay triangulation using the source coor¬ 
dinates. We used code developed by |BemaI| ( |1988| l which provides 
the coordinates of the vertices of the unique set of triangles over 
the FIRST area. We then computed the area covered by FIRST 
by adding up the areas of the individual triangles. The resulting 
area is 9,032.27 square degrees. Finally, for each FIRST source, 
we queried the SDSS database to see if the position of the sources 
was included in the survey. The result of the query was that 89.38% 
of FIRST sources fall in the SDSS-DR7 imaging area, implying an 
overlap of 8,073.04 square degrees. 













































































Radio quasars at z rsj /J_ 3 


Table 1. Steps in the pre-selection process 


3 SELECTION OF QSO CANDIDATES 


Pre-selection criterion 

Selected sources 

Starting FIRST-SDSSS sample (sj 1".5) 

222,517 

15 s: r s: 20.2 

74,853 

Point-like 

13,956 

Without mag error > 0.2 in all bands 

13,934 

Exclude sources with ’fatal’ error flags 

13,287 

Exclude sources with ’CHILD’ flag 

9,139 


2.2 Pre - selection criteria 


We matched each FIRST source, not flagged as possible sidelobe or 
nearby bright source (~ 3.6% of the sources in the catalog have this 
warning flag), with the closest optical object in the 'PhotoPrimary’ 
view of the SDSS DR7 catalogue within a l"5 radius. This radius is 
the same as used by C08, and is a compromise between complete¬ 
ness and efficiency. The adopted value is lower than the 2 , /0 radius 
used by SDSS in their algorithm for QSO selection. However, the 
excellent astrometry of FIRST and SDSS means that the peak in the 


distribution of optical/radio offsets occurs at about 072 (Schneider 
et al. 2010, Fig. 6), supporting our adoption of a l'/5 radius. From 


this match we obtained a starting sample of 222,517 sources. In 
this sample there is no selection by radio flux density or radio mor¬ 
phology other than the requirement that the radio source have at 
least a weak core component. The FIRST catalogue itself intro¬ 
duces several minor selection effects: the FIRST sensitivity limit is 
somewhat non-uniform over the sky, with small variations due to 
the observing strategy and large variations due to decreasing sensi¬ 
tivity in the vicinity of bright sources. However, the fraction of the 
survey area affected by sensitivity variations is small, less than 15% 
{ [Becker, White & Helfand| 1995fr . Another effect is that the FIRST 
survey limit of ^ 1.0 mjy refers to the peak flux density of sources 
rather than to the integrated flux density; consequently, extended 
sources with total fluxes greater than 1 mjy may not appear in the 
catalogue because their peaks fall below the detection threshold. 

From these 222,517 matches, we first selected the 13,956 star- 
like objects with 15.0 ^ r ^ 20.2, where r refers to SDSS 
psfMag_r, corrected for Galactic extinction according to|Schlegel,| 
|Finkbeiner & Dav!s| ( [1998| l. 

We then filtered the sample on the basis of several SDSS 
quality-control parameters used by others (e.g. ( [Richards et al.| 
|2002[ >) when selecting QSO targets for spectroscopy. Specifically, 
we rejected all objects with magnitude errors > 0.2 mag in all five 
bands, and any for which the SDSS ’fatal" error flags ’BRIGHT", 
"SATURATED", "EDGE" or "BLENDED" were set, indicating un¬ 
reliable photometry. This left 13,287 objects. 

Finally, in contrast with |Richards et aI7| < |2002) >, we rejected 
all objects with the "CHILD" flag set (another 4,148 sources), 
indicating objects obtained by de-blending an image flagged 
"BLENDED". This criterion ensures that only one optical object 
is associated with each radio source l |Carballo et al.[2006) . and we 
adopt it for consistency with C08. In this way we avoid introducing 
differences in the pre-selection that may change the final efficiency 
of the neural-network algorithm. This is the main source of incom¬ 
pleteness in our sample, as will be discussed in Section 8. 

This pre-selection process, summarised in Table[I] left us with 
9,139 star-like objects coinciding with FIRST radio sources. 


3.1 Neural Network algorithm 

The machine-learning technique used in this work to select our list 
of QSO candidates is described in earlier papers (C06, C08), and 
here we give only a brief summary. 

We used a supervised Artificial Neural Network (NN) algo¬ 
rithm of feed-forward type, suitable for solving classification prob¬ 
lems, and programmed using the Matlab Neural Network Toolbox 
software ®. A supervised NN is trained with samples of known 
classification, in order to learn how to distinguish between the 
classes. Only after the training has been carried out, can the trained 
NN be used to classify a new problem sample. In our case the 
classification task was formulated as a binary problem, the two 
classes being: (a) the target class, i.e. QSOs in the redshift range 
3.6 ^ z ^ 4.4, and (b) the "non-target" class, i.e. all other types of 
object. 

In a feed-forword NN, each input variable corresponds to a 
node in the so-called ‘input layer", Each of the input nodes has a 
weighted connection to every node in the next layer, called "the 
hidden layer". A node in the hidden layer forms a weighted sum 
of its inputs, and then passes the information to a second hidden 
layer that performs a similar processing. The number of hidden 
layers, like the other parameters of the NN, has to be optimised 
for the problem that the NN aims to solve. The weighted informa¬ 
tion passes through the layers of the NN to the last ’output’ layer, 
which performs a simple sum of its inputs, giving the output. 

In particular our NN was composed of the input layer, just one 
hidden layer, and an output layer y. The output y for the 7th object, 
with values in the range (0,1), is given by the non-linear function: 


V = 


1 + e~ 


( 1 ) 


with a 1 = wo + i Wjx'j, where (xr,X 2 ,-—,Xd) 1 are the in¬ 
put variables for object i. wo and (wi, W 2 , ...,Wd ), called bias 
and weights respectively, are the parameters fitted during the train¬ 
ing. This NN model is known as logistic linear discriminant. The 
adopted error function was the variance of the outputs: 


1 m 

e =- yy - r ) 2 (2) 

where m is the number of objects used for the training and T 
is the target value, set to 1 for the class of high-redshift QSOs and 
0 for the remaining sources, during the training. The optimal pa¬ 
rameters for the net. i.e. those minimising the error, were obtained 
using the Levenberg-Marquardt algorithm. This is a simple but ro¬ 
bust function for optimisation and it appears to be the fastest for the 
training of moderate-sized NNs ( jHagan & Menhaj| 1994} . 

The set of input variables adopted for the NN is the best set ob¬ 
tained in C08, i.e. a combination of optical magnitudes and colours 
(r,u — g,g — r,r — i,i — z) and radio-optical separation. The input 
variables were pre-processed, normalising their values to the range 
(—1, 1). No outliers required trimming. For this step the whole pre¬ 
selected sample was used, regardless of whether the source was 
spectroscopically classified and thus suitable as a training object, or 
not. In fact, the input variables of the new objects (i.e. the problem 
objects) presented to the trained net are expected to be normalised 
in the same way as the ones used in the training process. 
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3.2 Spectroscopic classification of the pre-selected sample 

Of the sample of 9,139 sources passing initial selection (Table 1), 
6,091 have spectra in the SpecObj view of SDSS-DR7. 5,348 of the 
latter are included in the 5th edition of the SDSS Quasar Catalogue 
(DR7 QSO Catalogue; |Schneider et al.|2010j >, which is also based 
on SDSS-DR7, but uses more stringent criteria for the classification 
of the objects as QSOs, in order to exclude dubious cases. 71 of the 
QSOs in this catalogue have redshifts in the range 3.6 ^ a ^ 4.4. 

The remaining 743 sources with spectra in DR7-SpecObj but 
not included in the DR7 QSO Catalogue are classified by SDSS 
as stars, galaxies, quasars and sources of ’unknown’ type. For all 
of these sources, a search was made in NED (the NASA/IPAC Ex- 
tragalactic Database), and none of them was classified there as a 
z ^ 3.6 QSO. We also visually inspected the DR7-SpecObj spectra 
of these sources to check if any of them could be a QSO in the red- 
shift range of interest here, but none of the objects, which include 
those with ’unknown’ spectra, had spectral features consistent with 
a high-z QSO. 

3,048 of the 9139 sources in the sample lack spectroscopic 
classification in DR7-SpecObj. Elowever, the DR7 QSO catalogue 
was compiled by inspecting all the SDSS spectra, not just those of 
the quasar candidates, and identified 115 of these as quasars (one 
of them being a high-z QSO in the redshift range of interest for 
our work) despite them not being automatically identified as such. 
Another 17 sources were identified as QSOs by C08, comprising 8 
high-z QSOs and 9 QSOs with redshift below 3.6. The remaining 
2,916 sources were checked in NED, and none of them had been 
spectroscopically classified as of March 2012. 

A total of 6,223 (= 6091 + 115 + 17) sources thus have a re¬ 
liable spectroscopic classification. Of these, 80 (= 71 + 1 + 8) are 
QSOs in the redshift range 3.6 ^ 4.4. These sources form 

the training sample, i.e. the sample used to train the neural network 
to distinguish high-z QSOs in our redshift range of interest from 
other objects. Quasars with 3.6 ^ ^ 4.4 play the role of the tar¬ 

get class.The non-target sources include stars, galaxies, QSOs with 
other redshifts and objects with spectra classed as ’unknown’ but 
lacking the features expected for our target sources. 

The 2916 sources without available spectra from SDSS-DR7 
or from the literature (as of March 2012), form the sample from 
which new QSO candidates in our redshift range of interest are 
selected using the trained NN. 


3.3 Training and testing of the NN 

Having defined the training sample, with 6223 spectroscopically 
classified sources, of which 80 are high-redshift QSOs, we are 
ready to train the NN and to test its performance as a classifier. 

The classification algorithm fitted by the NN provides for each 
source an output 0 ^ y ^ 1. The extreme values of 1 and 0 cor¬ 
respond respectively to sources with input variables more similar 
or less similar to those of the high-z QSO class. Objects with mea¬ 
sured y greater than some threshold value y c are candidate high-z 
QSOs. 

The performance of the trained NN can be expressed in terms 
of two basic parameters: efficiency (or reliability) and complete¬ 
ness. In our case the efficiency is the fraction of candidate high-z 
QSOs selected by the NN which are true high-z QSOs. The com¬ 
pleteness is the fraction of true high-z QSOs with y ^ y c , i.e., the 
fraction of true high-z QSOs selected as such by the NN. 

The performance of the NN is ideally tested with a sample 
of objects not used during the training. In our case, since the target 
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Figure 1 . Efficiency versus completeness, measured for the test sample. 
Each symbol corresponds to a given value of the threshold y c , ranging from 
0.9 to 0, in increments of 0.1, from left to right. The symbol with the error 
bars correspond to the adopted threshold y c = 0A. 

sample has only 80 objects, the ‘leave one out" method was applied, 
using all but one of the objects for the training, and the remaining 
one for the test. In total 6223 NNs were run, each of them provid¬ 
ing the ouput for a test object. These output values, for a sample of 
6223 test objects, were used to compute efficiency and complete¬ 
ness as a function of y c , and the results are shown in Figure[l] Since 
our purpose is to build a sample appropriate for statistical analysis, 
priority is given to completeness, accepting lower y values at the 
cost of lower efficiency. Choosing y c = 0.1 our NN classifier has 
an efficiency of 60% ± 9 and a completeness of 97% ±11 (errors 
assume Poisson statistics). 

3.4 High-redshift QSO candidates 

The set of 6223 trained NNs used for the testing was applied to 
the sample of 2916 sources without spectra, in order to find candi¬ 
date high-z QSOs. For each source we adopted the median of the 
6223 output values, and we selected as high-redshift QSO candi¬ 
dates those sources with j/ me d > 0.1. In this way 15 QSO candi¬ 
dates in the range 3.6 ^ ^ 4.4 were found (Table 2), out of an 

original set of 2,916 sources lacking spectroscopic identification. 


4 CHECKS OF THE QSO SELECTION 

Below we present spectroscopy of the 15 NN-selected high-z-QSO 
candidates (Section 4.1) and use SDSS DR9 to check the complete¬ 
ness of our selection (Section 4.2). In addition, we use our new data 
to assess the efficiency of high-z-QSO selection by C08 (Section 
4.3) and to check the completeness of the spectroscopic identifica¬ 
tion of high-z-QSOs in SDSS DR7 (Section 4.4). 

4.1 Spectroscopy of the 15 NN-selected candidates 

At the time, none of the 15 QSO candidates selected by our NN 
had spectra available in the literature. Therefore long-slit spectra 
of all 15 were obtained with the 2.5-m Nordic Optical Telescope 
(NOT) during the night of 25th March 2012, using ALFOSC (An- 
dalucla Faint Object Spectrograph and Camera) with grism #4, 
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which provides a dispersion of 3A per pixel. The spectral cover¬ 
age was 4000 — 9000 A and the resolution was 15A. The exposure 
time was about 1000s per source, delivering signal-to-noise ratios 
> 8 per pixel. The seeing was typically better than l’\3 FWHM. A 
spectrophotometric standard star was observed to correct for spec¬ 
tral response. After each target spectrum, an exposure of an arc- 
lamp was taken for wavelength calibration. Data were reduced us¬ 
ing standard IRAlQroutines. 

All 15 candidates were confirmed as QSOs, 7 of them (Fig. [2} 
in the desired range of redshift, 5 in the nearby range 3.1 ^ z < 3.6 
(Fig. 0- and the remaining 3 at lower redshifts (Fig.[4]l. For the last 
3 quasars at lower redshift, the NN-misclassification is caused most 
likely by confusion of Mg II emission with Lya emission. Table 
2 lists the coordinates, r magnitudes, y m ed values and redshifts 
for the 15 candidates. A note in column 6 indicates whether there 
is a spectrum of the object in the latest SDSS release DR9 (see 
next subsection). We report also the discovery of five new broad 
absorption line (BAL) QSOs (see Table 2). 

The observed efficiency of this selection is therefore ~ 47% 
(7 z > 3.6 QSOs out of 15 candidates). 

4.2 Check of completeness from SDSS DR9 

As discussed in earlier sections, our methodology is applied to 
a sample of sources detected in both SDSS and FIRST using a 
match radius of l'/5, and it makes use of photometric data in SDSS 
and of the radio-optical separation. Our pre-selection and selection 
methods are therefore based on the same variables as the SDSS 
target-selection algorithms. Therefore we checked if the more re¬ 
cent SDSS-DR9 spectroscopic catalogue provides spectra of any of 
the 2,916 sources lacking a spectral classification after SDSS-DR7 
was released. 

The ’SpecObj’ view of SDSS-DR9 was used for this pur¬ 
pose, giving the following results. 4 of the 15 candidates 
have spectra in DR9 (see Table [2J with redshifts very sim¬ 
ilar to those reported in this work. Of the remaining 2,901 
sources rejected by the NN as high-z QSO candidates, 451 have 
spectra in DR9 and 4 of them are classified as QSOs with 
3.6 ^ z ^ 4.4. However, examination of the spectra re¬ 
veals that all four are actually lower-redshift objects or stars. 
These objects are J075757.87+095607.56, J 101403.75+451053.27, 
J112742.74+363429.5, and J222758.13+003705.45. 

The fact that none of these 451 objects are classified in DR9 
as high-z QSOs is consistent with the estimated high completeness 
of our NN algorithm: ~ 97%. 

4.3 Assessment of the sample of high-z candidates in C08 

The selection of high-z QSOs in the present paper is based on more 
recent SDSS data releases than used by C08. However the FIRST- 
SDSS pairs were obtained in this work in same way as by C08, with 
similar criteria for the magnitude limits, the maximum radio-optical 
separation, optical morphology, and photometric quality. Also the 
classification procedure, aimed at the identification of new high-z 
QSOs, uses similar NN architecture, input parameter-set, NN train¬ 
ing and output parameter. 

1 IRAF is distributed by the National Optical Astronomy Observatories, 
which is operated by Association of Universities for Research in Astron¬ 
omy, Inc., under cooperative agreement with the National Science Founda¬ 
tion. 



Figure 5. Schematic overview of a check of the efficiency of our selection 
via comparison with DR9 spectroscopic information. See Section 4.2. The 
2916 objects (Section 3.2) without spectra in DR7, or in the literature, were 
classified by our NN. with 15 being selected by the NN as high-redshift 
candidates. (*) = targets from BOSS (Baryon Oscillations Spectroscopic 
Survey). 

In C08 we used a training sample comprising 52 QSOs with 
z ^ 3.6, all from the DR5 spectroscopic catalogue, and selected 58 
new candidates. In that paper 24 of the candidates were confirmed 
as z ^ 3.6 QSOs (17 from observations obtained for the paper 
and 7 from the literature or from the DR6 spectroscopic catalogue). 
16 sources were classified as other types of object, on the basis 
of observations by C08, or spectra from the literature or DR6. 18 
high-z QSO candidates remained unclassified at C08. 

11 of the 18 C08 candidates overlap with the 15 candidate- 
QSOs identified in this work. As result of our observations, these 
11 sources were classified as 5 high-z QSOs, 5 QSOs with 3.2 ^ 
z ^ 3.5 and a QSO at z = 1.40 (see Table[2]l. The remaining seven 
candidates of C08 consist of three QSOs included in the DR7 QSO 
Catalogue and in our training sample (SDSS 110946.44+190257.6 
with z=3.67, SDSS 123128.22+184714.3 with z=3.33 and SDSS 
124323.16+235842.2 with z=3.49) and four sources not selected 
by our classifier as being high-z QSOs. The spectra of three of 
these four sources were obtained in our observing programme at 
the NOT, yielding a QSOs at z = 3.44 and z = 1.4 and a late-type 
star (see Table[3] and figures[3]and[4]for spectra) . 

The spectroscopic observations of the 58 high-z candidates at 
C08 are now almost complete (57 out of 58). They yield an effi¬ 
ciency of 52 ± 9 per cent (30/58), highlighting the value of simple 
neural networks for this classification task. This efficiency is in rea¬ 
sonable agreement with the value obtained with the training sam¬ 
ple, i.e. with the expectation we had from the objects with known 
classification, 62 ± 9 per cent (C08). In addition, we note that a 
large fraction of the contaminants, 15 out of 28, are QSOs with 
3.1 ^ z < 3.6, close to the redshift threshold we adopted. 

4.4 Spectroscopic completeness of SDSS for high-z QSOs 

Several studies of the SDSS selection of QSOs ^ Richards et al. 
|2002||Croom et al.|2004||Richards et al.|2006||McGreer, Helfand & 

|White|2009[l suggest an overall completeness above 90%. The com- 
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Figure 2. NOT/ALFOSC spectra of NN-selected high-z candidates: QSOs with 3.6 ^ z ^ 4.4 
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Table 2. Sample of 15 FIRST-SDSS DR7 high-z QSO candidates selected by our NN 


RA 

DEC 

rAB 

Si.4GHz 

Vmed 

redshift 

DR9 

C08 

Notes 

(J2000) 


(mJy) 







(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

08:15:55.02 

+46:53:21.4 

19.89 

2.97 

0.12 

3.20 


yes 

BAL 

08:33:16.91 

+29:22:28.0 

20.13 

12.63 

0.32 

3.30 

yes 

yes 


08:57:24.33 

+ 11:05:49.2 

19.81 

1.91 

0.98 

3.71 




09:09:53.85 

+47:49:43.2 

19.90 

373.29 

0.22 

3.64 


yes 


09:14:36.23 

+50:38:48.5 

20.19 

47.98 

0.15 

3.62 


yes 


09:26:40.29 

-02:30:41.5 

19.82 

1.9 

0.12 

3.76 

yes 



10:29:40.93 

+ 10:04:10.9 

19.47 

2.81 

0.22 

3.40 


yes 


10:34:20.43 

+41:49:37.5 

20.12 

2.17 

0.30 

4.00 

yes 

yes 


11:33:00.71 

-04:11:58.5 

19.96 

9.64 

0.12 

3.39 



BAL 

11:51:07.42 

+50:15:58.6 

20.09 

1.69 

0.32 

3.40 


yes 

BAL 

12:05:31.73 

+29:01:49.2 

20.17 

1.51 

0.50 

3.44 


yes 

BAL 

12:13:29.43 

-03:27:25.7 

19.64 

23.37 

0.76 

3.67 


yes 


12:28:19.97 

+47:40:30.4 

19.32 

2.24 

0.46 

1.40 


yes 


12:44:43.07 

+06:09:34.6 

19.78 

1.29 

0.21 

3.76 

yes 

yes 


15:43:36.59 

+ 16:56:21.8 

18.97 

10.85 

0.12 

1.40 



LoBAL 


The columns give the following: (1) SDSS J2000 coordinates; (2) SDSS dereddened PSF r magnitude; (3) FIRST peak radio flux density; (4) NN output; (5) 
QSO redshift determined in this work; (6) indicates if the source has a spectrum in SDSS-DR9; (7) indicates if the source was previously selected as a high-z 
candidate by C08; (8) BAL - broad-absorption-line QSO; LoBAL - low-ionization broad-absorption-line QSO. 


Table 3. Spectra of three high-z QSO candidates from C08, not selected by our NN, but observed to complete the C08 sample 


RA 

DEC 

rAB 

Si.4GHz 

redshift 

DR9 

(J2000) 


(mJy) 



(1) 

(2) 

(3) 

(4) 

(5) 

08:48:18.88 

+39:38:06.0 

20.15 

1.28 

1.4 


10:58:07.47 

+03:30:59.6 

19.91 

4.18 

3.44 

yes 

12:04:07.84 

+48:45:48.0 

19.97 

3.96 

M-star 



The columns give the following: (1) SDSS J2000 coordinates; (2) SDSS dereddened PSF r magnitude; (3) FIRST peak radio flux density; (4) QSO redshift 

or spectral classification; (5) indicates if the source has a spectrum in SDSS-DR9. 


pleteness in general decreases with increasing redshift and decreas¬ 
ing brightness, and it is particularly inefficient for 2.2 < 2 < 3 
where quasar and star colours are very similar. In particular, |Mc-| 
|Greer, Helfand & White| i |2009fr studied the completeness of quasar 
selection at redshift z > 3.5 and magnitude i < 20.2, and found 
for the SDSS target algorithm a completeness of « 86%, in good 
agreement with the 85% derived in|Richards et al.|lj2006fr. 


The analysis in this paper and in C08 identifies 15 QSOs with 
3.6 ^ z ^ 4.4, missed by SDSS-DR7 (two have spectra in SDSS- 
DR9, which uses BOSS). This allows us to estimate the incomplete¬ 
ness of the SDSS-DR7 selection, after allowing for the fact that the 
spectroscopic area of the SDSS survey is ~ 95% of the imaging 
area. 4 of our 15 QSOs 3.6 ^ z ^ 4.4 lie outside the SDSS spec¬ 
troscopic plates, so the estimated incompleteness in SDSS is 11 
QSOs out of 83 (87 minus 4) in this redshift range, i.e. ~ 13% ± 4, 
in good agreement with estimates by McGreer, Helfand & White 
( [2009] ) and |Richards et al.| l |2002| l. 


In Fig.[6]we plot the g — i,r — z colour-colour diagram of the 
entire training sample together with the 15 3.6 ^ z ^ 4.4 QSOs 
missed by SDSS and identified by us. Splitting the sample between 
QSOs with 3.6 ^ z ^ 4.4 and all other sources, we note that 
our QSO-selection method is sensitive even at the boundary of the 
two samples, demonstrating the effectiveness of learning-machine 
techniques, when compared with simple colour-cut criteria. 



g-i 

Figure 6. Colour-colour diagram of the 6233 sources with spectra in DR7, 
with green crosses representing QSOs with 3.6 ^ z ^ 4.4. and red crosses 
representing other objects. The new 15 QSOs with 3.6 ^ z ^ 4.4 (8 from 
C08, 7 from this work, all missed by SDSS) are plotted as black diamonds. 
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Figure 3. NOT/ALFOSC spectra of NN-selected high-z candidates: QSOs with 3.0 ^ z < 3.6 


5 FINAL SAMPLE 

In Table[4]we present the final sample of 87 QSOs with 3.6 ^ 2 ^ 
4.4 satisfying our selection criteria. This sample includes 72 QSOs 
(see Section 3.2) from the 5th quasar catalogue (Schneider et al.| 
|2010| > plus 15 QSOs revealed by our neural-network (8 from C08 
and 7 from this work, see Section 4.1). The magnitide limit for our 


sample is tab = 20.2. To convert to absolute magnitudes M r , a 
/C-correction (Humason, Mayall & Sandage| 1956} is required. 


Following the convention of |Hogg et al.| (2002} , the K- 
correction between a bandpass R used to observe a source at red- 
shift 2 and the same bandpass in the rest-frame, is: 
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Figure 4. NOT/ALFOSC spectra of NN-selected high-z candidates: QSOs with z < 3 


m R = M r + D m {z) + K r (z) 


(3) 


where Dm(z) is the distance modulus calculated from the lumi¬ 
nosity distance Dl as Dm = 5logia rriR is the apparent 

magnitude, and Mr is the absolute magnitude. 

An accurate Tf-correction, including the contribution of emis¬ 
sion lines, can be computed by convolving a typical QSO spectrum 
at different redshifts with the filter response (e.g. |Cristiani & Vio| 
1 1990 and Wisotzki||2000| the latter based on optical/UV spectra 
i Elvis et al.|1994). Following this approach, we calculated the 


from 


K-correction by convolving the Vanden Berk et al. (2001 1 compos¬ 
ite quasar spectrum with the SDSS r-filter: 


K = 2.51og 10 


fo°° F(\)S(\)dX ' 

/ 0 °°F((A)/(l + *))S(A)dA. 


(4) 


where F( A) is the measured intensity per unit wavelength, 
and S( A) is the r-band filter response. The resulting K-correction 
is shown in Fig. [7] and in Table[5| 

To convert from M r to the commonly used monochromatic 
absolute AB magnitude at 1450A, M 1450 , we assume the canonical 
power-law spectral energy distribution with spectral index a„ = 
— 0.5 and hence: 


AT 1450 = M r + 2.5a„log 1 


/1450 A 


= M r + 0.791 (5) 


310 \6231 °A t 

where 623lA is the effective wavelength of the SDSS r filter. We 
use this spectral index rather than that derived by |Vanden Berk et al.| 
<|200l) ( a v = —0.44) to allow a direct comparison with other au¬ 
thors. We al so avoid using the composite spectrum in Vanden Berk| 


et al. 


(2001) to convert from 6231A to 1450A rest-frame flux den¬ 


sity since above H/3 there is a significant contribution from stellar 



Figure 7. Computed A'-correction in r. A tabulation is available on-line, 
and a sample is shown in Table 5. 


light from low-redshift quasar hosts, which we effectively eliminate 
by extrapolating the power-law to red wavelengths. 


6 COMPLETENESS 

Below we estimate the completeness (number of QSOs selected di¬ 
vided by the actual number of QSOs in this range of redshift) of 
the sample of 87 radio-loud QSOs with 3.6 ^ z ^ 4.4 . This 
sample will be used (in Section 7) to calculate the luminosity func¬ 
tion. There are several sources of incompleteness: exclusion of op¬ 
tical images of poor quality (Section 6.1); radio-survey incomplete¬ 
ness and missed radio-optical identifications (Section 6.2); and in- 
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Table 4. Final sample of 87 radio-loud QSOs with 3.6 ^ z 4.4 (continued on next page) 


NAME 

RA 

DEC 

rAB 

<7 f- 

Redshift 

Si.4GHz 

M r 

logi 0 Pi.4GHz 

ID 


(J2000) 




(mJy) 


(W/Hz) 


(1) 


(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

J015339.61-001104.9 

01:53:39.61 

-00:11:05.0 

18.83 

0.022 

4.19 

4.82 

-28.59 

26.42 

SDSS 

J030025.23+003224.2 

03:00:25.23 

+00:32:24.2 

19.68 

0.025 

4.18 

7.75 

-27.74 

26.62 

SDSS 

J072518.26+370518.3 

07:25:18.27 

+37:05:18.4 

19.60 

0.020 

4.33 

26.72 

-27.94 

27.19 

WHT 

J074711.14+273903.3 

07:47:11.15 

+27:39:03.4 

18.35 

0.025 

4.15 

1.08 

-29.04 

25.76 

SDSS 

J074738.49+133747.3 

07:47:38.49 

+13:37:47.3 

19.35 

0.015 

4.17 

7.18 

-28.05 

26.59 

WHT 

J075113.04+312038.0 

07:51:13.05 

+31:20:38.0 

19.73 

0.020 

3.76 

5.84 

-27.29 

26.42 

SDSS 

J075122.35+452334.2 

07:51:22.36 

+45:23:34.2 

20.20 

0.033 

3.61 

1.18 

-26.73 

25.69 

SDSS 

J080710.74+131739.4 

08:07:10.74 

+13:17:39.4 

20.00 

0.026 

3.73 

48.20 

-27.00 

27.32 

SDSS 

J081009.95+384757.0 

08:10:09.95 

+38:47:57.1 

19.62 

0.018 

3.94 

26.68 

-27.57 

27.11 

SDSS 

J082323.32+155206.8 

08:23:23.32 

+15:52:06.8 

19.30 

0.018 

3.78 

74.93 

-27.74 

27.53 

SDSS 

J083322.50+095941.2 

08:33:22.50 

+09:59:41.2 

18.69 

0.016 

3.73 

122.52 

-28.31 

27.73 

SDSS 

J083808.46+534809.8 

08:38:08.46 

+53:48:09.9 

19.94 

0.032 

3.61 

8.47 

-27.00 

26.54 

SDSS 

J083946.22+511202.8 

08:39:46.22 

+51:12:02.9 

19.31 

0.016 

4.39 

40.50 

-28.28 

27.38 

SDSS 

J084044.19+341101.6 

08:40:44.19 

+34:11:01.6 

19.78 

0.020 

3.89 

13.64 

-27.35 

26.81 

SDSS 

J085257.12+243103.1 

08:52:57.12 

+24:31:03.2 

19.46 

0.016 

3.62 

157.30 

-27.48 

27.81 

SDSS 

J085501.81+182437.7 

08:55:01.82 

+18:24:37.7 

19.96 

0.020 

3.96 

9.38 

-27.25 

26.66 

SDSS 

J085724.32+110549.1 

08:57:24.33 

+11:05:49.2 

19.81 

0.017 

3.71 

1.91 

-27.17 

25.92 

NOT 

J085944.06+212511.1 

08:59:44.06 

+21:25:11.2 

18.74 

0.015 

3.70 

23.54 

-28.24 

27.01 

SDSS 

J090254.16+413506.5 

09:02:54.17 

+41:35:06.5 

20.12 

0.023 

3.69 

1.12 

-26.85 

25.68 

WHT 

J090953.85+474943.2 

09:09:53.85 

+47:49:43.2 

19.90 

0.020 

3.64 

373.29 

-27.05 

28.19 

NOT 

J091436.23+503848.5 

09:14:36.23 

+50:38:48.5 

20.19 

0.028 

3.62 

47.98 

-26.75 

27.30 

NOT 

J091824.38+063653.3 

09:18:24.38 

+06:36:53.4 

19.76 

0.022 

4.19 

25.87 

-27.66 

27.15 

SDSS 

J092640.28—023041.4 

09:26:40.29 

-02:30:41.5 

19.82 

0.021 

3.76 

1.90 

-27.20 

25.93 

NOT 

J092832.87+184824.3 

09:28:32.88 

+18:48:24.4 

17.54 

0.016 

3.77 

8.67 

-29.49 

26.59 

SDSS 

J093714.48+082858.5 

09:37:14.49 

+08:28:58.5 

18.58 

0.015 

3.70 

3.47 

-28.40 

26.18 

SDSS 

J094003.03+511602.7 

09:40:03.03 

+51:16:02.7 

19.00 

0.014 

3.60 

12.90 

-27.94 

26.72 

SDSS 

J100012.26+102151.8 

10:00:12.26 

+10:21:51.9 

19.54 

0.023 

3.64 

21.87 

-27.41 

26.96 

SDSS 

J101747.75+342737.8 

10:17:47.76 

+34:27:37.9 

20.00 

0.031 

3.69 

2.83 

-26.98 

26.08 

SDSS 

J103055.95+432037.7 

10:30:55.95 

+43:20:37.7 

19.84 

0.025 

3.70 

31.29 

-27.14 

27.13 

SDSS 

J103420.43+414937.5 

10:34:20.43 

+41:49:37.5 

20.12 

0.029 

4.00 

2.17 

-27.12 

26.03 

NOT 

J103446.54+110214.4 

10:34:46.54 

+11:02:14.5 

18.81 

0.025 

4.27 

1.15 

-28.68 

25.81 

SDSS 

J103717.72+182303.0 

10:37:17.72 

+18:23:03.1 

19.81 

0.023 

4.05 

13.59 

-27.48 

26.84 

SDSS 

J105121.36+612038.0 

10:51:21.37 

+61:20:38.1 

18.92 

0.020 

3.69 

6.90 

-28.06 

26.47 

SDSS 

J105756.25+455553.0 

10:57:56.26 

+45:55:53.1 

17.45 

0.022 

4.14 

1.10 

-29.93 

25.77 

SDSS 

J110147.88+001039.4 

11:01:47.89 

+00:10:39.4 

20.18 

0.026 

3.69 

192.10 

-26.79 

27.92 

SDSS 

11105:43.86+255343.1 

11:05:43.87 

+25:53:43.1 

20.09 

0.026 

3.75 

1.69 

-26.92 

25.87 

SDSS 

11109:46.44+190257.6 

11:09:46.44 

+19:02:57.6 

20.05 

0.024 

3.67 

7.22 

-26.91 

26.49 

SDSS 

11110:55.21+430510.0 

11:10:55.22 

+43:05:10.1 

18.58 

0.024 

3.82 

1.21 

-28.50 

25.74 

SDSS 

11117:01.89+131115.4 

11:17:01.90 

+13:11:15.4 

18.28 

0.018 

3.62 

28.00 

-28.66 

27.07 

SDSS 

11117:36.33+445655.6 

11:17:36.33 

+44:56:55.7 

20.05 

0.026 

3.85 

24.33 

-27.05 

27.05 

SDSS 

11123:39.59+291710.7 

11:23:39.60 

+29:17:10.8 

19.46 

0.017 

3.77 

3.68 

-27.57 

26.22 

SDSS 

11125:30.49+575722.7 

11:25:30.50 

+57:57:22.7 

19.43 

0.036 

3.68 

2.52 

-27.54 

26.03 

SDSS 

11127:49.45+051140.5 

11:27:49.45 

+05:11:40.6 

19.14 

0.012 

3.71 

2.34 

-27.84 

26.01 

SDSS 

11129:38.73+131232.2 

11:29:38.73 

+13:12:32.3 

18.78 

0.026 

3.61 

1.77 

-28.16 

25.86 

SDSS 

11133:30.91+380638.1 

11:33:30.91 

+38:06:38.2 

19.71 

0.025 

3.63 

1.39 

-27.23 

25.76 

SDSS 


The columns give the following: (1) SDSS object-ID (2) SDSS J2000 coordinates; (3) SDSS dereddened PSF r magnitude; (4) error in PSF r magnitude as 
given in SDSS; (5) QSO redshift determined in this work or from SDSS; (6) FIRST peak radio flux density; (7) absolute r magnitude ; (8) radio luminosity at 
1.4 GHz; (9) indicates the source of the data from which the redshift was first obtained (the two WHT sources are from|Benn et al.|2002|. 


completeness of selection by the neural-network algorithm (Section 
6.3). 

6.1 Optical image quality 

Due to the sensitivity of our NN to data of poor photometric qual¬ 
ity, we discarded at the pre-selection stage (Section 2.2) objects 
having ’fatal’ error flags or magnitude errors larger than 0.2 in all 
five bands (Richards et al. 2002) or flagged CHILD. 

Incompleteness due to the exclusion of fatal and non-fatal 
photometric errors ( |Richards et al.|2002| define as ’non-fatal’ er¬ 
rors: some empirical combination of SDSS flags generally asso¬ 


ciated with poor de-blends of complex objects) during the SDSS 
selection of QSOs-candidates, was discussed by [Richards et al.| 
( |2006j ) who applied a global 5% correction. In previous evaluations 
of this selection effect, (Croom et al. 2004 ) suggest 6% incomplete¬ 
ness for objects with 17.5 i 18.5, and [Van den Ber k et al.| 
l |2005j ) estimated an incompleteness of 3.8% for point-like objects 
with i < 19.1 . 

To quantify the incompleteness of our selection due to the ex¬ 
clusion of CHILD objects, we evaluated the fraction of such objects 
amongst 15 ^ r ^ 20.2 QSOs in the 5th SDSS Quasar Catalogue, 
which did not exclude CHILD objects. The fraction that we derive 
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Table 4. (continued) 


NAME 

RA 

DEC 

TAB 

(X 7 - 

Redshift 

Sl.4GHz 

M r 

logioPl.4GHz 

ID 


(J2000) 




(mJy) 


(W/Hz) 


(1) 


(2) 

(3) 

(4) 

(5) 

(6) 

(?) 

(8) 

(9) 

J113729.42+375224.2 

11:37:29.43 

+37:52:24.2 

20.18 

0.032 

4.17 

2.70 

-27.22 

26.16 

SDSS 

J115045.61+424001.1 

11:50:45.61 

+42:40:01.1 

19.88 

0.019 

3.87 

1.82 

-27.23 

25.93 

SDSS 

J115449.36+180204.3 

11:54:49.36 

+18:02:04.4 

19.61 

0.024 

3.69 

37.99 

-27.36 

27.21 

SDSS 

J120447.15+330938.7 

12:04:47.15 

+33:09:38.8 

19.23 

0.025 

3.62 

1.10 

-27.71 

25.66 

SDSS 

J121329.42—032725.7 

12:13:29.43 

-03:27:25.7 

19.64 

0.027 

3.67 

23.37 

-27.33 

27.00 

NOT 

J122027.95+261903.5 

12:20:27.96 

+26:19:03.6 

18.13 

0.017 

3.70 

34.34 

-28.85 

27.17 

SDSS 

J123142.17+381658.9 

12:31:42.17 

+38:16:58.9 

20.18 

0.030 

4.14 

20.44 

-27.19 

27.03 

SDSS 

J124054.91+543652.2 

12:40:54.91 

+54:36:52.2 

19.74 

0.023 

3.94 

14.89 

-27.44 

26.86 

SDSS 

J124209.81+372005.6 

12:42:09.81 

+37:20:05.6 

19.34 

0.018 

3.84 

644.79 

-27.75 

28.47 

SDSS 

J124443.06+060934.6 

12:44:43.07 

+06:09:34.6 

19.78 

0.024 

3.76 

1.29 

-27.24 

25.76 

NOT 

J124658.82+120854.7 

12:46:58.83 

+12:08:54.7 

20.01 

0.024 

3.80 

1.07 

-27.05 

25.69 

SDSS 

J124943.67+152707.0 

12:49:43.67 

+15:27:07.1 

19.34 

0.019 

3.99 

1.75 

-27.90 

25.94 

SDSS 

J130348.94+002010.5 

13:03:48.94 

+00:20:10.5 

18.89 

0.019 

3.65 

1.08 

-28.06 

25.66 

SDSS 

J130738.83+150752.0 

13:07:38.83 

+15:07:52.1 

19.72 

0.027 

4.11 

3.44 

-27.63 

26.26 

SDSS 

J131242.86+084105.0 

13:12:42.86 

+08:41:05.0 

18.52 

0.014 

3.74 

4.41 

-28.49 

26.29 

SDSS 

J131536.57+485629.0 

13:15:36.58 

+48:56:29.1 

19.76 

0.025 

3.62 

9.94 

-27.18 

26.61 

SDSS 

J132512.49+112329.7 

13:25:12.49 

+11:23:29.8 

19.32 

0.022 

4.41 

69.39 

-28.28 

27.61 

SDSS 

J134854.37+171149.6 

13:48:54.37 

+17:11:49.6 

19.13 

0.021 

3.62 

2.10 

-27.81 

25.94 

SDSS 

J135406.89—020603.2 

13:54:06.90 

-02:06:03.2 

19.17 

0.018 

3.72 

709.05 

-27.82 

28.49 

SDSS 

J135554.56+450421.0 

13:55:54.56 

+45:04:21.1 

19.36 

0.021 

4.09 

1.48 

-27.98 

25.89 

SDSS 

J140635.66+622543.3 

14:06:35.67 

+62:25:43.4 

19.73 

0.020 

3.89 

11.03 

-27.41 

26.72 

WHT 

J140850.91+020522.7 

14:08:50.91 

+02:05:22.7 

19.08 

0.017 

4.01 

1.47 

-28.17 

25.87 

SDSS 

J142209.70+465932.5 

14:22:09.70 

+46:59:32.5 

19.72 

0.022 

3.81 

10.56 

-27.35 

26.68 

SDSS 

J142326.48+391226.2 

14:23:26.48 

+39:12:26.3 

20.15 

0.024 

3.92 

6.63 

-27.01 

26.50 

SDSS 

J143413.05+162852.7 

14:34:13.06 

+16:28:52.7 

19.86 

0.022 

4.19 

4.90 

-27.56 

26.42 

SDSS 

J143533.78+543559.2 

14:35:33.78 

+54:35:59.2 

20.04 

0.025 

3.81 

93.26 

-27.02 

27.63 

SDSS 

J144542.76+490248.9 

14:45:42.76 

+49:02:48.9 

17.32 

0.009 

3.87 

2.51 

-29.80 

26.07 

SDSS 

J144643.36+602714.3 

14:46:43.37 

+60:27:14.4 

19.79 

0.033 

3.78 

1.87 

-27.25 

25.92 

SDSS 

J145329.01+481724.9 

14:53:29.01 

+48:17:24.9 

20.12 

0.030 

3.77 

4.42 

-26.91 

26.30 

WHT 

J150328.88+041949.0 

15:03:28.89 

+04:19:49.0 

17.96 

0.017 

3.66 

122.70 

-29.00 

27.72 

SDSS 

J150643.80+533134.4 

15:06:43.81 

+53:31:34.5 

18.94 

0.022 

3.79 

13.97 

-28.11 

26.80 

SDSS 

J151146.99+252424.3 

15:11:46.99 

+25:24:24.3 

19.95 

0.024 

3.72 

1.24 

-27.04 

25.73 

SDSS 

J152028.14+183556.1 

15:20:28.14 

+18:35:56.2 

19.82 

0.021 

4.12 

6.26 

-27.54 

26.52 

SDSS 

J153336.13+054356.5 

15:33:36.14 

+05:43:56.5 

19.84 

0.020 

3.94 

27.55 

-27.34 

27.13 

SDSS 

J161716.49+250208.1 

16:17:16.49 

+25:02:08.2 

19.98 

0.023 

3.94 

1.01 

-27.20 

25.69 

SDSS 

J 161933.65+302115.0 

16:19:33.65 

+30:21:15.1 

19.53 

0.025 

3.81 

3.88 

-27.53 

26.25 

SDSS 

J163708.29+091424.5 

16:37:08.30 

+09:14:24.6 

19.57 

0.018 

3.75 

9.51 

-27.45 

26.62 

WHT 

J163950.52+434003.6 

16:39:50.52 

+43:40:03.7 

17.95 

0.017 

3.99 

23.78 

-29.28 

27.07 

SDSS 

J164326.24+410343.4 

16:43:26.24 

+41:03:43.4 

20.10 

0.025 

3.87 

63.40 

-27.02 

27.47 

SDSS 

J222814.40-085525.7 

22:28:14.40 

-08:55:25.7 

20.18 

0.034 

3.64 

1.99 

-26.77 

25.92 

WHT 

J223535.59+003602.0 

22:35:35.59 

+00:36:02.1 

20.15 

0.027 

3.87 

5.06 

-26.97 

26.38 

NOT 

J235022.40-095144.3 

23:50:22.40 

-09:51:44.4 

19.67 

0.021 

3.70 

6.37 

-27.31 

26.44 

WHT 


The columns give the following: (1) SDSS object-ID (2) SDSS J2000 coordinates; (3) SDSS dereddened PSF r magnitude; (4) error in PSF r magnitude as 
given in SDSS; (5) QSO redshift determined in this work or from SDSS; (6) FIRST peak radio flux density; (7) absolute r magnitude ; (8) radio luminosity at 
1.4 GFlz; (9) indicates where the source redshift was first obtained (the two NED sources are from|Benn et al.|2002}. 


z em K correction 

0.01 . 0.0075 

0.02. 0.0151 

0.03 . 0.0241 

0.04. 0.0412 

0.05 . 0.0720 

0.06 . 0.1236 

0.07. 0.1619 

0.08 . 0.1851 

0.09. 0.1959 

0.10. 0.1999 

0.11. 0.1985 

0.12. 0.1939 


Table 5. K -correction in the SDSS r band (Fig. 7). A portion of the table 
is shown here. The full table is available on-line 


in this way is 27%. The net completeness due to exclusion of these 
two types of object is therefore 69% (0.95 x 0.73). 

These fractions are in approximate agreement with the statis¬ 
tics of Table 1, which indicate 4.8% incompleteness due to ’fatal’ 
errors (669 rejected out of a total of 13956 sources), and 31% in¬ 
completeness due to the exclusion of sources flagged as CHILD 
(4148 out of 13287), i.e. a net completeness of 66% (0.95 x 0.69). 

For the analysis here, we adopt an intermediate estimated 
completeness of 68%. 

We do not apply any correction for objects misclassified in 
SDSS as having galaxy morphology but being star-like objects. 
Tests of random samples of sources with available spectra indicate 
that this source of incompleteness is < 0.03%. 
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Figure 8. Cumulative distribution of FIRST flux densities for those QSOs 
in our sample of 87 with 3.6 ^ z ^ 4.4, and with radio flux density < 10 
mJy. 


6.2 Radio incompleteness 

Two sources of incompleteness arise from radio selection using the 
FIRST survey. 


6.2.1 Incompleteness of FIRST sur\>ey 

The completeness of the FIRST survey as a function of flux den¬ 
sity has been estimated for SDSS quasars and is given in fig. 1 of 
Jiang et al. l |2007} (and is discussed elsewhere, e.g. by jPrandoni] 
et al.pOOT) . At integrated FIRST flux density > 3 mJy FIRST 
is > 96% complete and the completeness declines with decreasing 
flux density. From the results of these studies we assume the follow¬ 
ing completeness function, q(S), where S is the FIRST integrated 
flux density in mJy: 


0.50, 


1.25 

0.75, 

1.25 

<5^2 

0.85, 

2 < 

S ^ 3 

0.95, 

3 < 

5^5 

1, 

S > 

5 


In Fig. [8] we show for our sample the cumulative distribution 
of integrated radio flux densities, for objects fainter than 10 mJy. 

Applying q(S) to our final sample we obtain a completeness 
of 84 per cent (87 QSOs detected and ~ 16 missed), for the flux 
density limit of 1 mJy. 


6.2.2 Match radius 

To obtain our radio-optical sample we sought simple one-to-one 
matches within radius F’,5. We used this criteria for consistency 
with C08, who adopted this on the basis that more than 99% of 
FIRST-APM QSOs with 3.8 sj z ^ 4.5, E sC 18.8 and S 1 AG Hz > 
1 mJy fall within this radius |Vigotti et al~~| ( |2003^ . 

However, a simple one-to-one match between FIRST and 
SDSS will miss double-lobe QSOs without detected radio cores. 


De Vries, Becker & White (2006) found that for a sample of 5, 515 
FIRST-SDSS QSOs with radio morphological information within 
450”, the fraction of FIRST-SDSS double-lobe QSOs with unde¬ 
tected cores is 3.7%. Since the starting samples of SDSS QSOs 
in < [de Vries, Becker & White[20()6j and in this work obey similar 
SDSS selection criteria, we used this value to correct for this source 
of incompleteness. 


6.3 Incompleteness of selection by the neural-network 

The completeness of our neural network classifier was estimated as 
97 per cent, from the testing on known high-z QSOs presented in 
Section 3.3. The classifier selected 15 high-z QSO candidates and 
rejected 2901. Of the later, 451 have now spectroscopy from DR9 
and none of them was identified as a high-z QSO, confirming this 
high level of completeness (see Section 4.2). 


6.4 Net completeness 

The net completeness of our sample of 87 RL-QSOs at 3.6 ^ a ^ 

4.4 (Section 5), is the product of each of the completeness terms 
discussed above: acceptance of only those candidates with high- 
quality photometry (Section 6.1,68% completeness); completeness 
of the FIRST radio survey (Section 6.2.1, 84% for our sample); 
acceptance of only optical-radio matches within a given radius and 
exclusion of extended sources (Section 6.2.2, 99% and 96.3%); and 
the completeness of the NN selection algorithm (Section 6.3,97%). 

Multiplying these four terms together, the net completeness 
for our sample is 53%. 


7 BINNED LUMINOSITY FUNCTION 
7.1 Method 


Using the final sample of 87 QSOs listed in Section 5 and correct¬ 
ing for incompleteness as discussed in Section 6, we compute the 
binned Quasar Luminosity Function (QLF) in the redshift range 
3.6 sj z 4.4. 


The binned QLF is usually calculated using the classical 
XjVmax method (Schmidt 1968 Maccacaro et al. 1991|;[Ellis et al. 


1996 1 , or its generalized version (usually known as ^ V a ) ap¬ 


plied to samples comprising subsamples with different flux limits 


( jAvni & Bahcall|l980| l. 

The V,nax method is an unbiased |Felten|1976| l non-parametric 
estimator of the space density. It is commonly used to fit models of 
the LF. since it has the advantage that it does not assume any un¬ 
derlying model. Even when the model LF is fitted to the unbinned 
data (for example in the maximum likelihood technique of |Marshall| 
jet al.|1983) it is often used before performing the fit to observe the 
overall behavior of the LF. 

However |Page & Carreraj (2000) demonstrated that the 
"^2 1/Vrnax estimator introduces significant errors for objects close 
to the flux limits of the survey. An alternative method proposed by 
|Page & Carreraj l [2000j > is superior and partially corrects for this 
source of error, although implicitly assumes a uniform distribution 
of the sources within each bin ( jCroom et al.|2 009 ; Miyaji, Hasinger 
& Schmidt 2001). The variation of the LF within a bin can be par¬ 
ticularly critical at the steep bright end of the QSO LF. Instead, we 
used a modified version of the jPage & Carreraj ( |2000| l method that 
does not make use of the uniform-distribution assumption and is 
still model-independent. 
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To illustrate the method used in this paper, we start with a 
brief overview of the l/Vmax and the Page & Carrera l |2000| 
methods. The luminosity function is defined as the number of 
objects per unit of comoving volume, per unit of luminosity. A 
naive approach to the calculation of space density in an interval 
[Li,L 2 ] x [Vi(zi), V 2 (z 2 )] of luminosity and redshift, centered 
upon values L* and z*, would be to simply count the number of 
objects N within the interval considered: 


& Carrera ( 2000 ) approach, (b) in the classical 1/Vma.x case, and 


(c) for the methodology used here. 

The statistical uncertainty 5$ is calculated for each bin i as 


»• - m <12) 

Where No is the actual number of objects in the bin and 
5No = y/No (Poisson statistics). The formula is easily derived, 
when the space density <f>i is assumed to be approximately: 


VKl (7) 

The 1/Vmax method, first proposed by!Schmidt (1968j) takes into 
account the fact that in flux-limited samples there is a higher prob¬ 
ability to observe a bright source than a faint one. Thus, the count 
of sources N is replaced with a sum of probabilities: 


N Np-f 
V e L e V e L e 


(13) 


where ( V E L e ) is an equivalent space-luminosity area, and N 
is the corrected number of quasars in bin i. The error in the com¬ 
pleteness factor / is assumed to be « 0 . 


E 

i 


Vo 

Vm ax,i 


( 8 ) 


where Vo is the volume over which we are computing the lumi¬ 
nosity function,and V m ax,i is the maximum volume at which the 
source could be observed and still be included in our sample. In 
this way, the computation of the LF becomes 


HL,z) 


1 N l 

1 E 1 

i 


A L 


AV r 


max,i 


(9) 


|Page & Carreraj ( |2000| ) noted that the limit in apparent magnitude 
of the survey bounds the region of integration. In particular, for 
a given bin in redshift, Li and L 2 should be replaced by the ac¬ 
tual luminosity limits (L m i n , L max ) as determined by the intersec¬ 
tions with the limiting-magnitude curves of the survey. Therefore, 
< 1 >(L, z) is calculated as 


HL,z) 



N 

r x max(L) ( dv/dz)dz dL 

J z min 


( 10 ) 


where z m i n is the bottom of the redshift interval and z m ax{L) is 
the highest possible redshift for an object of luminosity L within 
the considered bin A z. This approach takes into account the real 
area of integration but implicitly assumes uniform distribution of 
sources over the bin. 

In order to minimise this bias, we calculate the maximum ac¬ 
tual integration area determined as in [Page & Carrera| l |2000j > but for 
each source in the bin, so as not to lose the V max information for 
individual sources. Then we sum over all the sources in the bin. In 
this way we do not count the number of sources over an area larger 
than that of the actual survey and at same time, within a single bin, 
we weight sources by luminosity. Finally, 4?(L, z) is calculated as 


S(L,*) = E 


£W W ,<( d V/dz)dzdL 


(ID 


The difference between this methodology and that of |Page &| 
|Carrera| ( |2000| > is negligible in the case of a large sample of QSOs 
uniformly distributed in (M,z) space, but becomes critical for small 
samples not distributed uniformly in each bin, as is the case here. In 
Fig ,|9]we show the volume-luminosity space available to an object 
in a bin intersected by a limiting-magnitude curve, (a) in the |Page| 


7.2 The QSO Luminosity Function 

The QLF was computed using two bins in redshift, 3.6 - 4.015 and 
4.015 - 4.415, and 11 bins in optical absolute magnitude starting 
with M r = —26.6 and with AM = 0.3. In Fig.[l0]we plot M r vs. 
redshift for our sample of QSOs; the dotted grid shows the bins in 
magnitude and redshift used to compute the QLF. The curves show 
the upper and lower limiting apparent magnitude r of our selection. 
The top and side panels show the marginal distributions in redshift 
and absolute magnitude, respectively. 

Since our complete sample results from two surveys with dif¬ 
ferent flux limits, the maximum redshift at which a source can be 
observed may be different in the two surveys. The most efficient 
way to combine areas with different flux limits is to assume that 
each object could be found in any of the survey areas for which it is 
brighter than the corresponding flux limit. This is ‘coherent’ addi¬ 
tion of samples, in the language of |Avni & B ahcall (1980) . There¬ 
fore, since our survey is a radio-optical survey, for each source we 
chose the smaller of optical and Z m ax,radio- 

For each bin we applied the completeness corrections as ex¬ 
plained in Section 6 and computed the weighted number of QSOs 
in the bin. Tableland Figure[TT]show the resulting QLF for radio- 
loud QSOs. In Fig. m we plot separately the luminosity function 
for the two bins of redshift, at z ~ 3.8, and z ~ 4.2. We also show 
the best-fit slopes which will be discussed in Sect. 8.3. 


8 DISCUSSION 

Below we use the binned LF calculated from the previous section 
to derive the space density of RL QSOs at 3.6 ^ z ^ 4.4 (Section 
9.1). We then derive the space density of the entire population (i.e 
radio-loud + radio-quiet) of QSOs at this redshift (Section 8.2) us¬ 
ing two independent estimates of the radio-loud-fraction (Sections 
8.2.1, 8.2.2). Finally we derive the slope of the LF of RL QSOs at 
a ~ 3.8 and at z ~ 4.2, and compare our results with those of 
other authors (Section 8.3). 

Our sample includes QSOs with optical luminosities M r < 
—26.6 4=> M 1450 < —25.8 (see equation [3^ We adopt the 
definition of radio-loudness used by Gregg et al.] ( |1996| ), i.e. 
log .Pi. 4 ,GHz > 25.5. All the sources included in our sample (i.e. 
<Si.4GHz > 1 mjy) meet this criterion and are therefore radio-loud. 
Due to the limit in flux density, the minimum radio luminosity of 
sources included in our sample depends on redshift. In particular, 
using a radio spectral index of a r = —0.3, at the lowest redshift of 
our sample, i.e. 2 = 3.6, the FIRST flux-density limit corresponds 
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Figure 9. Volume-luminosity space used to calculate the contribution to the LF from a single object (the black dot) in a given bin intersected by the line 
L = Ln rn (z), i.e. the minimum detectable luminosity of an object at redshift 2 . The available space (grey shaded area) is shown for a binned LF calculated 
using: (a) the Page & Carrera implementation , (b) the classical 1/Vmax method, (c) our methodology. 
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Figure 10. Luminosity-redshift diagram for the complete sample of 87 radio loud QSOs. The dotted lines show the limits of the bins ( ALAz ) used to compute 
the luminosity function. Solid lines represent the upper (r = 15) and lower (r = 20.2) limits of the survey. Upper and left panels show marginal histograms 
of redshift and absolute magnitude, respectively (see Section 7.2) 


to radio luminosity log Pia,gb.z(W/ Hz) > 25.61. For z — 4, it 
corresponds to log Pia,ghz(W/Hz) > 25.7 and for z = 4.4 to 
log Pi.a,gbm{W/Hz) > 25.77. 

8.1 The space density of RL QSOs at 3.6 < z < 4.4 

Starting with the binned luminosity function determined in Sec¬ 
tion 8.2 we calculate the space density of RL QSOs with 


optical luminosity M 1450 < —25.8 and radio luminosity 

log 1 ^ 1 . 4 , GHz (IV I Hz) > 25.5, in two shells of redshift. The first 
shell has median z ~ 3.8(3 .6 ^ z ^ 4.015) and the second shell 
has median a « 4.2 (4.015 ^ 2 ^ 4.415). Integrating the binned 
LF (Table 6 ), the space densities of QSOs are therefore: 
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Table 6. Binned luminosity function for FIRST-SDSS quasars at redshift ~ 4 


Redshift 

(1) 

M r 

(2) 

<Wxio- 9 ) 

(3) 

o>(xl0 9 ) 
(4) 

RLF 

(5) 

</>(x 10 -8 ) 
16) 

till 

(7) 

Nq 

(8) 

NQ cor r 

(9) 

N r 

(10) 

3.8 

-26.75 

3.129 

1.400 

3.17(±1.60)% 

9.870 

1 

5 

9.43 

0 

3.8 

-27.05 

3.778 

0.890 

3.58(±1.82)% 

10.559 

1 

18 

33.96 

0 

3.8 

-27.35 

3.242 

0.744 

4.04(±2.07)% 

8.025 

0 

19 

35.85 

1 

3.8 

-27.65 

1.240 

0.469 

4.56(±2.36)% 

2.717 

0 

7 

13.21 

1 

3.8 

-27.95 

1.260 

0.476 

5.15(±2.69)% 

2.447 

0 

7 

13.21 

1 

3.8 

-28.25 

0.849 

0.379 

5.80(±3.01)% 

1.463 

0 

5 

9.43 

0 

3.8 

-28.55 

0.679 

0.339 

6.53(±3.50)% 

1.040 

0 

4 

7.55 

0 

3.8 

-28.85 

0.170 

0.170 

7.35(±4.00)% 

0.231 

0 

1 

1.89 

0 

3.8 

-29.15 

0.340 

0.240 

8.26(±4.58)% 

0.411 

0 

2 

3.77 

0 

3.8 

-29.45 

0.170 

0.170 

9.26(±5.24)% 

0.183 

0 

1 

1.89 

0 

3.8 

-29.75 

0.170 

0.170 

10.39(±6.01)% 

0.163 

0 

1 

1.89 

0 

4.2 

-27.05 

1.920 

1.920 

3.04(±1.60)% 

6.316 

1 

1 

1.89 

0 

4.2 

-27.35 

0.657 

0.464 

3.44(±1.81)% 

1.910 

1 

2 

3.77 

0 

4.2 

-27.65 

0.917 

0.410 

3.89(±2.07)% 

2.358 

0 

5 

9.43 

0 

4.2 

-27.95 

0.550 

0.318 

4.39(±2.36)% 

1.254 

0 

3 

5.66 

0 

4.2 

-28.25 

0.367 

0.260 

4.95(±2.69)% 

0.741 

0 

2 

3.77 

0 

4.2 

-28.55 

0.367 

0.260 

5.58(±3.07)% 

0.658 

0 

2 

3.77 

0 

4.2 

-29.15 

0.183 

0.183 

7.07(±4.01)% 

0.260 

0 

1 

1.89 

0 

4.2 

-30.05 

0.183 

0.183 

10.02(±6.03)% 

0.183 

0 

1 

1.89 

0 


The columns give the following: (1) median redshift of the bin, (2) median absolute magnitude M r of the bin, (3) space de nsity l (Mpc~ 3 
mag -1 (x 10 -9 )) of the radio-loud QSOs , (4) error on the space density, cr^Y. 10 -9 ), (5) radio-loud fraction calculated using Jiang et al. |2007 , (6)) 
space density </> ( Mpc ~ 3 mag~ 1 (x 10 —8 )) of the QSOs (RQ+RL), (7) indication if the bin is intersected by the limiting magnitude curve (1 = yes and 0 = 
not), (8) actual number of QSOs in the bin, (9) the corrected number of QSOs in the bin after applying the completeness corrections, (10) Number of QSOs 

in the bin limited by the radio-flux limit. 



Figure 11. The luminosity function derived from our sample of 87 radio- 
loud QSOs with log 10 -Pi. 4 ,GHz(W/Hz) > 25.7 . The luminosity func¬ 
tion for z 3.8 is shown with black points, while red squares show the 
luminosity function for z ~ 4.2. Black and red lines are the best fit-slopes. 


p{z « 3.8, A/1450 < -25.8)rl = 4.51 ± 0.61 Gpc -3 
p{z « 4.2, Mi 450 < — 25.8)rl = 1.54 ± 0.63 Gpc -3 

From a sample of radio QSOs obtained by cross-matching the 
FIRST radio survey and the Automatic Plate Measuring Facility 
catalogue of POSS I, | Vigotti et al.| ( j2003) measured the space den¬ 
sity at 3.8 ^ z ^ 4.5 of optically-luminous (A/1450 < —26.9) 
radio-loud QSOs and obtained p(z ~ 4.1, A/1450 < — 26.9)rl = 
0.99 ± 0.28 Gpc -3 . We recalculated the space density and opti¬ 
cal luminosities using our adopted cosmology (noted at the end 
of Section 1), obtaining p(z m 4.1, A/1450 < — 27.1)rl = 


0.66 ± 0.18 Gpc 3 . By integrating our binned LF in the interval 
A/1450 < —27.0, we obtain: 


p{z « 3.8, A/1450 < — 27.0)rl = 1.09 ± 0.24 Gpc -3 
p{z « 4.2, M1450 < — 27.0)rl = 0.50 ± 0.16 Gpc -3 


in good agreement withjVigotti et al.l j2003) (see Fig.|12[l and con¬ 
sistent with a linear decrease of space density with increasing red¬ 
shift. 

Using a sample of QSOs obtained by cross-matching FIRST 
and SDSS-DR6, |McGreer, Helfand & White] <2009) calculated a 
binned luminosity function in the redshift range 3.5 ^ a ^ 4.0. 
These authors used the same starting surveys as we did and a sim¬ 
ilar range of redshift, but they calculated the LF only for QSOs 
with radio-loudness parameter R > 70. The R parameter is an¬ 
other common criterion for distinguishing between radio-quiet and 
radio-loud AGN. It is defined ( [Kellermann et al.|1989HStocke et al.| 
|1992| ) as the rest-frame ratio of the monochromatic 6-cm (5 GHz) 
and 2500A flux densities. Generally, objects are considered to be 
RL for R > 10. 

The space density calculated in |McGreer, Helfand & White| 
( 20091 for A/1450 < —26.1 is p(z = 3.75, A/1450 < 

— 26 . 1 )r> 7 o = 1.38 ± 0.59 Gpc -3 . The cosmology parameters 
used by |McGreer, Helfand & White| <2009) are the same that we 
use. 


At redshift ss 4 our definition of radio-loudness is very close 
to the common definition R > 10, but we needed to re-calculate 
the LF using a subsample of RL-QSOs with R > 70 in order to 
compare our LF with |McGreer, Helfand & White[ ( |2009) . To calcu¬ 
late the //-parameter for our sample of QSOs, we used a v = —0.5 
(in agreement with |McGreer, Helfand & White 2009 1 to transform 
the flux from 5'i.4GHz to S'sghz- We follow |Oke & Gunn] <1983) 
when converting from magnitude to luminosity (2500 A). In this 









































16 D. Tuccillo et al. 



Figure 12. The black points show the integrated luminosity function for RL 
QSOs with M 1450 < —27.0. The points for z ~ 3.8 and z ~ 4.2 were 
obtained in this work, showing good agreement with the point obtained for 
z ~ 4.1 in Vigotti et. (2003). For comparison, blue points show the space 
density for the entire population of QSOs and for M 1450 < —27.6, as 
found in Richards et al. (2006). See Section 8.1. 


way, we obtain 



Figure 13. Cumulative luminosity functions at z = 3.8 (black circles) and 
z = 4.2 (red circles). Squares represent densities derived by | Vigotti et al. j 
(2003) (blue square) and Carballo et al. (2006 1 (purple square) at redshifts 
z ~ 4.1 and z ~ 4, respectively. 


We therefore derive below the space density of all QSOs in 
two different ways: assuming constant RLF (Section 8.2.1); and us¬ 
ing a recently-determined redshift- and luminosity- dependent RLF 
(Jiang et al,|2007) (Section 8.2.2). 


p{z = 3 . 8 ,Mi 45 o < —26.1 )r> 70 = 2.49 ± 0.36 Gpc -3 (14) 


which is a factor 1. 8 (2a ) higher than the value p = 1.38 ± 
0.59 Gpc -3 found by 


McGreer, Helfand & White 


(2009). This dif¬ 


ference may in part be ascribed to the higher completeness of our 
NN selection, and in part to the smaller FIRST-SDSS matching ra¬ 
dius used by [McGreer, Helfand & WhTte|(2009| . which will exclude 
some quasars. 

In Fig. 13 we show the cumulative luminosity functions for 
the two redshift bins (i.e. z = 3.8 and z = 4.2). Each point of the 
cumulative function is the space density p(< M 1450 ) as a function 
of absolute magnitude. The two functions can be compared with 
previous results by |Vigotti et al,|(2003) and |Carballo et al.] ( |2006| >, 
at redshifts z ~ 4.1 and z ~ 4, respectively. As expected, due to 
the evolution of space density with redshift, these last two values 
lie between our determinations. 


8.2.1 For constant RLF 

From our binned luminosity function we derive the space density 
for A /1450 < —26.1 (equivalent to the limit of A/; < —27.6 in 
|Richards et al.|2006| > as: 


p(z = 3.8, M1450 < — 26.1)rl = 3.57 ± 0.44 Gpc -3 
p(z = 4.2, M1450 < — 26.1)rl = 1.54 ± 0.63 Gpc -3 

In |Vigotti et al.|p003) > the RLF is assumed almost constant and is 
estimated as: 

RLF(Mi 4 so < -26.9) = 13.3% 

Therefore, for the total QSO population we obtain: 


8.2 Total space density of QSOs at 3.6 ^ z ^ 4.4 


From the space density of the radio-loud QSO population we can 
roughly test the predictions of Radio-Loud-Fraction (RLF) over 
this range of redshift, by comparing the implied total space den¬ 
sities with measurements of space density from the literature. 

It has long been known that between 5% and 15% of all 
quasars are radio-loud (e.g., Kellermann et al 7| 1 1 989) |Urry &| 
|Padova ni 1995; Ivezic et al. 2002). However, some authors con¬ 
clude that the RLF does not change significantly with redshift (e.g., 
|Goldschmidt et al.ll999[|Stem et al.|2000||Cirasuolo et al.l2003| or 


luminosity (e.g., Bischof & Becker|1997[|Stem et al.|2000|Vigotti| 


|et al.|2003 ), while others find that the RLF decreases with ir 


ing redshift (e.g.. Peacock, Miller & Longair|1986 Miller, 

^eacock] 

& Mead 19901 Visnovsky et al. 1992 

Schneider et al. 1992 

and de- 

creasing opical luminosity (Padovani 

1993, La Franca et a 

. 1994); 


or that it evolves non-monotonically with redshift and luminosity 
(e.g. |Hooper et al.|1995| l. 


p(z = 3.8, A/1450 < -26.1) = 26.8 ± 3.3 Gpc -3 
p(z = 4.2, A/1450 < -26.1) = 11.6 ±4.7 Gpc -3 

in good agreement with the space densities derived from Richards 
et al. (2006): 

p(z = 3.75, Mi 450 < -26.1) = 29.0 ± 2.0 Gpc -3 
p(z = 4.2, Mi 450 < -26.1) = 13.6 ± 1.3 Gpc -3 


8.2.2 For redshift- and luminosity-dependent RLF 

Jiang et al. (2007) l use a sample of more than 30,000 optically se¬ 
lected QSOs from the SDSS to study the evolution of the RLF as 
a function of redshift and luminosity. They find that the RLF of 
QSOs decreases with increasing redshift and decreasing luminos¬ 
ity, according to: 
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1° S (1 - RLF) =bo + b dog(l + z) + 6 m(A/ 2500 + 26) (15) 

where A/2500 is the absolute magnitude at rest frame 2500 
A. The parameters b 0 , b z and 6 m depend on the value of radio 
loudness and for R > 10 ; they are b 0 = —0.132 ± 0.116, b z = 
—2.052±0.261, 6m = —0.183±0.025. M2500 is calculated from 
M r as 

M 25 oo = M r + 2.5a, log (16) 

We used the Jiang et al. ( |2007|> formula to obtain for each bin 
the corresponding value of the RLF (column (5) in Table 5). The 
RLF lies in the range 3%-10%, and increases with decreasing 0 rl- 
Applying the corresponding RLF to each bin and integrating we 
obtain: 


p(z x 3.8, A /1450 < -26.1) = 81.7 ± 31.7 Gpc " 3 
p(z « 4.2, A/1450 < -26.1) = 41.0 ± 31.1 Gpc " 3 

This is a factor ~ 3 higher than the results from |Richards et ah] 
( |2006| >, but still within 2 cr, due to the large errors in the luminosity 
function and in the RLF. In particular, the errors on the RLF at this 
redshift and magnitude are ~ 50%. Given that our determination of 
the radio-loud luminosity function agrees reasonably well with |Mc-| 
| Greer, Flelfand & White) ( |2009| ) and with |Vigotti et al.| ( |2003| ), this 
discrepancy cannot be attributed solely to a possible overestimation 
of our luminosity function but may also be due to to a systematic 


underestimation of the RLF in Jiang et al. (2007 

. The large quoted 

errors invite caution when using the Jiang et a” 

(20071 formula to 


determine the fraction of radio-loud quasars at high redshifts. 


8.3 The bright-end slope of the luminosity function for RL 
QSOs 

The QLF is usually well fitted by a double-power-law parametri- 
sation that takes into account the redshift (e.g. |Pei|1995| Peterson 
[l997l|Boyle et al.|2000||Croom et al.|2004||Richards et al.|2006| l: 


<j>(L,z) 


_ 07 L* _ 

(L/L*(z))- + {L/L*(z))~P 


(17) 


where a, /3, 0* , L* are the faint-end slope, the bright-end slope, 
the normalisation of the luminosity function, and the characteris¬ 
tic break luminosity, respectively. This model, if a > P, can be 
approximated by 


0 OC 



ifL>L* 
if I, < L* 


(18) 


As already mentioned we calculated the LF in terms of opti¬ 
cal luminosity in two bins of redshift. We compare our results and 
the best-fit slope with those of Rich ards et al.| (|2 006|> for the en- 
tire population of QSOs. and with the results of McGreer, Helfand 
|& White| ( |2009) for RL QSO with R > 70. The limiting magni¬ 
tude of the QSOs samples used by these authors (A/1450 < —26.1 
for Richards et al.| ( |2006) and |McGreer, Helfand & White[ l |2009) l) 
was considered bright enough and far from the break luminosity to 
approximate the LF by a single power law oc L fi . 

This kind of approximation led in recent years to a long debate 


about an apparent flattening of the bright-end slope for z > 4, after 
it was noticed in early high-redshift surveys (jSchmidt, Schneider &| 
|Gunn|1995]|Fan et al.|2001) . These authors showed that the slope 
at z > 4 had a value /3 « —2.5, much shallower that the one seen 
at z < 2.2 (/? = —3.3, Croom et al. 2004). This flattening was 
then confirmed by |Richards et al. | ( [2006( 1 who used a large, homo¬ 
geneous QSO sample from the SDSS-DR3 extending to 2 = 5. 
At higher redshift the constraints are weaker as they come from 
small samples, but in general they do not confirm a continued flat¬ 
tening of the slope with increasing redshift. In fact, |Willott et al.| 
( |2010) >. combining the CFHQS (Canada-France High-z quasar sur¬ 
vey) with the more luminous SDSS sample, derived the QLF from 
a sample of 40 QSOs at redshifts 5.74 < a < 6.42 and found 
—3.8 < f3 < —2.3. At redshift z ~ 6,|Jiang et al.| (2008} find 
P = -3.1 ± 0.4 using QSOs from SDSS Stripe 82. 


Evolution of the shape of the QLF with redshift (changes in 
the slopes or in the location of the break luminosity) provides one 
of the fundamental observational constraints to the growth of super- 
massive black holes (SMBHs) over cosmic time. Assuming that 
brighter AGN have more-massive black holes, the flattening of the 
bright-end would be a remarkable indication of a downsizing of the 
SMBHs at high redshift. Downsizing was reported also by X-ray 
surve ys {U eda et a l.|2003[ [Basinger, Miyaji & Schmidt|2005| see 
also |Brusa et al.|2009| i! 

On the other hand, recent work by |Shen & Kelly| ( j2012[ l and 
|McGreer et al.| f2013 ) aims to fill the gap in the QLF between 
z ~ 3.5 and z ~ 6 , with the purpose of testing the flattening of the 
bright-end slope at z > 3. |Shen & Kelly|{2012ft constrain the lumi¬ 
nosity function by Bayesian modeling and using an homogeneous 
sample of SDSS-DR7 QSOs at z = 0.3 — 5. The results of|S hen &| 
|Kelly| < {20 1 2) and |Richards et al.| {2006) are, in general, in good 
agreement, finding that the curvature of the LF changes signifi¬ 
cantly beyond z = 3. However, |Shen & Kelly|{2012} suggest that 
the apparent flattening of the slope appears to be more related to a 
strong evolution of the break luminosity than a change in the bright- 
end slope. A similar conclusion is drawn by |McGreer et al.| ( |2013| >, 
who find no evidence for an evolution in the bright-end slope at 
A /1450 < —26 for a sample of QSOs with 4.7 ^ z ^ 5.1. On the 
other hand, |McGreer et al.| < [2013| > find evidence of strong evolu¬ 
tion in the break luminosity, as it brightens from A/i 45 o ~ —25.4 
at z = 2.5 to A /1450 « —27.2 at z = 5. They conclude that 
this evolution could flatten the bright-end slopes for surveys where 
the faint limit is near the break luminosity. |McGreer et al.| {2013) 
compared different models for the evolution of the QLF normal¬ 
ization and break luminosity. Eventually they found a good fit of 
their data with recent results from the literature, using a modified 
version of a luminosity evolution and density evolution (LEDE) 
model proposed by |Ross et al.|{20l2) . In particular the evolution 
of the break luminosity in this model is log-linear (up to z ~ 5), 
with a break luminosity that brightens with redshift. This modified 
LEDE model predicts that for z ~ 3.8 the break luminosity would 
be A/4450 ~ —26.2 and for z ~ 4.2 it would be A/i 450 ~ —26.4. 


If we approximate the LF by a single power law oc /7, we find 
that in the first bin of redshift, z ~ 3.8, our best-fit slope is P = 
—2.3 ± 0.2. As shown in Fig |14|our best-fit is in good agreement 
with the slope found in Richards et al. {2006)7i = — 2.4±0.1), and 
with that found by [McGreer, Helfand & White) < (2009| > (P = —2.2± 
0.2). For the RL-QLF calculated in our second bin of redshift, i.e. 
z ~ 4.2, we re-binned the LF using AM = 0.6, in order to reduce 
the statistical noise. In this way, as shown in Fig[l5]the best-fit slope 
is P = —2.0 ± 0.4. This result is consistent with the result found 
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by (Richards et al.j |2006| for the entire population of QSOs, i.e. 
/3 = -2.2 ±0.1. 

Our determinations of the bright-end slope for the RL popula¬ 
tion of QSOs at 2 ~ 3.8 and a ~ 4.2 are consistent with the flat¬ 
tening (between these redshifts-bins) of the bright-end slope found 
in |Richards et al.H2006) for z ^ 4, which [McGreer et al,| \2() 1 3) 
suggest is due to a bias resulting from a single-fit power law in a 
region near the break luminosity. 

As we have quasars with luminosities near or below the pre¬ 
dicted break luminosity, we repeat the fit but excluding those 
points. In the first bin at redshift z ~ 3.8, we exclude the two 
fainter points. In this way we obtain a slightly steeper best-fit slope 
and a larger error: /3 = —2.4 ± 0.3. This fit is shown in Fig| 1 4| as 
a dashed line. In the second bin of redshift, z ~ 4.2, we exclude 
the faintest point, obtaining again a small increase of the slope: 
ft = —2.1 ± 0.4. This fit is shown again as a dashed line in Fig|15| 

In light of the results from |McGreer et ak| ( |2013| l, our data do 
not strongly constrain the slope of the bright end nor the exact loca¬ 
tion of the break luminosity, especially considering the large errors 
of the brighter bins of the LF. Nevertheless this simple derivation is 
consistent with the results of |McGreer et al.| l j20 13} . 

In summary, our results are in good agreement with those of 
|Richards et al.|{2006| > and |McGreer, Helfand & White| ( |2009) . This 
result in itself is not trivial, because we are comparing different 
populations of QSOs in this range of redshift. In particular, we are 
comparing our results with the whole population of QSOs (by com¬ 
paring with |Richards et al.|2006[ > and with a population of RL QSOs 
where the radio-loudness is defined differently (being R* > 70 for 
the RL sample of |McGreer, Helfand & White|2009) . We therefore 
have indications of a certain homogeneity of the QLF regardless of 
the differences in radio-loudness. 

On the other hand, since we can’t constrain the bright and the 
faint-end slopes, and we don’t have an estimate of the break lumi¬ 
nosity, we can’t conclude that the noted consistency of the slopes 
imply also consistency of all the parameters. Any differences of 
the slopes, or a different value of the break luminosity, could point 
to different density evolution of the RL and the RQ populations. 
In fact, |Jiang et aL] ( j2007[ > express the dependence of the RLF on 
optical luminosity as « L 0 5 , implying Prlqso ~ /3qso + 0.5, 
which is consistent with the differences of slopes that we find be¬ 
tween the first (z ~ 3.8) and the second bin (z ~ 4.2) of redshift. 
Also, |Balokovic et al.| \20\2) and |Kratzer| (2014) ) find evidence 
that, at high redshift, the radio-loudness distribution of quasars is 
not a universal function, and likely depends on redshift and/or op¬ 
tical luminosity. We therefore need a sample of RL QSOs at fainter 
luminosities to constrain the faint-slope, and a larger survey area 
to extend the bright end of the luminosity function and thus, deter¬ 
mine the break luminosity. 


9 CONCLUSIONS 

We construct a sample of high-redshift radio-loud QSOs 3.6 ^ z ^ 
4, and use it to measure the luminosity function and space density 
of QSOs in this range of redshift. Our principal conclusions are: 

(i) We show (Section 3) that a simple neural network can be 
used to select high-redshift QSOs from radio-optical surveys, with 
97% completeness and 60% efficiency. 

(ii) With the aid of the neural-network, we construct a sample of 
87 radio-loud QSOs at redshift ~ 4. Of the various sources of in¬ 
completeness in the optical and radio surveys (Section 6 ), exclusion 
of SDSS ’CHILD’ images is the main cause of our incompleteness 



Figure 14. Black points show the luminosity function derived for z ~ 3.8, 
for RL QSOs with logioPi 4 qr^(W/H z) > 25.7 . For comparison, red 
points show the luminosity function calculated by McGreer et al. (2009), for 
RL QSOs with radio-loudness R > 70. and blue points the LF as calculated 
by Richards et al.(2006) for the entire population of QSOs. Best-fit slopes 
are /3 = —2.3 for our LF (black line) ; /3 = —2.2 for McGreer et al. 2009 
(red line), /3 = —2.4 for Richards et al.(2006) (blue line). For our LF we 
obtain /3 = —2.3 (black line) and /3 = —2.4 after excluding the two fainter 
points (dashed line). See Section 8.3. 



Figure IS. Black points show the luminosity function derived for z ~ 4.2, 
for RL QSO with lo(ju) P\ .4 .gI I AW /Hz) > 25.7 . For comparison, blue 
points show the LF as calculated by Richards et al. (2006) for the entire 
population of QSOs in the same bin of redshift. Best-fit slope is /3 = —2.2 
for Richards et al. (2006) (blue line). For our LF we obtain /3 = —2.0 
(black line) and /3 = —2.1 after excluding the faintest point (dashed line) 
See section 8.3 


relative to the SDSS selection of QSO candidates. But when ap¬ 
plied to non-CHILD objects, our neural-network algorithm detects 
~ 97% of the high-z QSOs, while SDSS only detect ~ 85% of 
them. 

(iii) We determine the optical luminosity function for radio-loud 
QSOs in two redshift bins, 3.6 ^ z < 4.0 and 4.0 ^ z ^ 4.4 (Sec¬ 
tion 7), and measure the total comoving density of QSOs in these 
two redshift ranges (Fig. 12), obtaining a result consistent with that 
of |Vigotti et al. | ( [2003) at luminosities M 1450 < —27.6. We also 
find good agreement between our cumulative luminosity functions 
(Fig. 13) and that measured by|Vigotti et al.| <[2003[> and|Carballo 
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|et al.| ( [2006^ , which determine the space density at intermediate 
redshifts. 

(iv) Assuming a radio-loud fraction of 13.3% jVigotti et al.| 
|2003f > we estimate the total comoving density of QSOs (Section 
8.2.1). The derived density of QSOs at z ~ 4 is consistent with 
that of [Richards et al.| < |2006| >. Alternatively (Section 8.2.2), using 
the redshift- and luminosity-dependent radio-loud fraction found 
by [Jiang et al. ( |2007| >, we measure a total comoving density of 
QSOs a factor 3 higher than measured by Richards et al. (20061. 
However, this result is significantly affected by the large error bars 
on the formula assumed for the radio-loud fraction (RLF). 

(v) We determine the slope of the luminosity function in two 
bins of redshift (Section 8.3). In the lower-redshift bin (z = 3.8) 
we found /3 = —2.3 ± 0.2, consistent with [Richards et ak| l |2006} 
and |McGreer, Helfand & White] |2009] l. In the higher redshift bin 
( z = 4.2) we find a slope /3 — 2.0 ± 0.4 consistent with [Richards] 
|et al.|p006| . Values of the slope consistent with our determination, 
have been interpreted as a flattening of the bright end slope for 
the high-z QSOs population, but has recently been re-interpreted 
as the result of a strong evolution of the break luminosity for high- 
z QSO |McG reer, Helfand & White| < |2009}.The consistency of our 
results with [Richards et al.| < |2006} > and lMcGreer, Helfand & White| 
l |2009} suggests a similar evolution for both radio-loud and radio¬ 
quiet populations. Our results can be also interpreted as sugges¬ 
tive of a flattening of the bright-end slope from z ~ 2 to z ~ 4, 
for the radio-loud population only. If confirmed, this implies an 
evolution of the density of super-massive black holes associated 
with radio-loud QSOs, in the sense that they were more abundant 
at z ~ 4. However, to clarify the evolution of the RL population 
relative to that of the whole population of QSOs, more observa¬ 
tional constraints are needed, especially at redshifts above 4. The 
candidate-selection approach described here is now being applied 
to FIRST-SDSS-UKIDSS surveys to search for QSOs at z > 4.5 
(Tuccillo, McMahon & Gonzalez-Serrano, in prep.) 
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